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Abstract. We study the fractional quantum Hall phases of a pseudospin-1/2 Bose 
gas in an artificial gauge field. In addition to an external magnetic field, the gauge 
field also mimics an intrinsic spin-orbit coupling of the Rashba type. While the spin 
degeneracy of the Landau levels is lifted by the spin-orbit coupling, the crossing of 
two Landau levels at certain coupling strengths gives rise to a new degeneracy. We 
therefore take into account two Landau levels, and perform exact diagonalization of 
the many-body Hamiltonian. We study and characterize the quantum Hall phases 
which occur in the vicinity of the degeneracy point. Notably, we describe the different 
states appearing at the Laughlin filling, v = 1/2. While for this filling incompressible 
phases disappear at the degeneracy point, denser systems at v = 3/2 and v = 2 are 
found to be clearly gapped. For filling factors v = 2/3 and v — 4/3, we discuss the 
connection of the exact ground state to the non-Abelian spin singlet states, obtained 
as the ground state of k + 1 body contact interactions. 



PACS numbers: 67.85.Fg,73.43.-f 



X 



Keywords: Quantum Hall states. Artificial gauge fields. Spin-orbit coupling. Ultracold 
bosons. 



FQH states of a Bose gas with spin-orbit coupling 



2 



1. Introduction 

Cold atoms provide versatile and highly controllable quantum systems. This makes 
them particularly useful as quantum simulators, with plenty of applications which stem 
from their complicated many-body behavior. Cold atoms have been used to explore 
important open problems in solid state physics. A noteworthy example is the physics 
emerging from strongly correlated electronic systems, which has fascinated researchers 
for decades. Due to its topological nature and the exotic anyonic excitations, fractional 
quantum Hall systems (FQH) are particularly interesting candidates for a quantum 
simulation with cold atoms [IJ EJ [3j E] . 

A major difference between electrons and atoms, however, is the atoms' lack of 
electric charge. This makes them insensitive to external magnetic fields, which would 
hinder exploring the vast field of quantum Hall physics. Different proposals to overcome 
this have been made: In principle, rotating the atomic system is sufficient to mimic the 
Lorentz force by mean of the Coriolis force [5j [6] and bring the system into quantum Hall 
phases [TJ, [8J [9]. However, experimental attempts to reach this regime by rotating the 
cloud have not been entirely successful [TQl HH H21 UHl HI] . Alternatively, the magnetic 
field can be mimicked by dressing the atoms with laser fields in such a way that the lowest 
dressed state obtains a space-dependent phase factor. Properly chosen, this phase factor 
acts as if the atoms were moving like charged particles in a magnetic field pJ3 [THJ [17] . A 
Raman dressing of 87 Rb atoms has been used to generate an artificial electric field [T8] . 
Later on, this scheme has successfully mimicked an artificial magnetic field and visibly 
induced vortices into the atomic cloud [19J. Artificial gauge fields have also been realized 
in optical lattices, by means of laser dressing [201 El] or shaking of the lattice |22j . 

These schemes allow for simulating the behavior of charged spinless (or spin- 
polarized) particles. The spin degree of freedom of electrons might play a role in 
fractional quantum Hall setups, as the strong magnetic field does not necessarily fully 
polarize the electron spin due to the small gyromagnetic ratio in many solid materials. 
Allowing the atoms to be in two internal states, one might try to mimic spin- 1/2 particles 
like electrons, but again the lack of charge will cause significant differences between this 
pseudospin, and the real spin of the electrons. The latter interacts with a magnetic 
field which is induced by the motion of the charged electron through the electric field of 
the nuclei. This effect, called Rashba spin-orbit coupling [23], is absent in cold atoms. 
Ways to overcome this have been proposed recently. The laser dressing of atoms can 
be designed in such a way that it not only implements the effect of a magnetic field, 
but also of a spin-orbit coupling [2U [25], [26] . Mathematically, spin-orbit coupling and 
magnetic fields can be treated on the same footing: While the magnetic field can be 
taken into account through a minimally coupled vector potential A, i.e. a term p ■ A 
in the Hamiltonian, the spin-orbit coupling has the form p ■ er, where er = (a x , o~ y , o~ z ) T , 
a vector of Pauli matrices. This makes the spin-orbit coupling equivalent to gauge 
fields, which in the case of spin- 1/2 belong to the SU(2) group, and are thus in general 
referred to as non-Abelian gauge fields. In pioneering experiments, spin-orbit coupling 
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has been synthesized experimentally for 87 Rb fZT\ [25], using a Raman dressing similar 
to the one of Refs. [HI [19]. These experiments have attracted a great deal of attention, 
and the impact of spin-orbit coupling on the condensation properties of a Bose gas has 
been studied theoretically [291 USB ED E21 [33j E31 [35] . Proposals to synthesize spin-orbit 
coupling for atoms in optical lattices [361 EZ] also motivate the study of Bose-Hubbard 
models with spin-orbit coupling [381 EH HO]- 

A different scenario is encountered in the presence of both a spin-orbit coupling 
and a strong artificial magnetic field. The latter induces strong correlations, producing 
fractional quantum Hall phases which in turn may be qualitatively affected by the spin- 
orbit coupling. Its influence could perfectly be revealed by a cold atom experiment 
which combines the artificial spin-orbit coupling and the artificial magnetic field, and 
allows to freely tune the coupling strength. The closest connection to the electronic 
fractional quantum Hall effect can be made by considering fermionic atoms with dipolar 
interactions [31] • 

In this paper, we consider the experimentally more feasible scenario of bosonic 
atoms with short-range interactions. They are confined to a plane which, assuming 
periodic boundary conditions, assumes a toroidal geometry. The bosons are subjected 
to an artificial magnetic field perpendicular to the plane, and we assume a Rashba 
coupling between an effective pseudospin-1/2 degree of freedom and the external motion 
of the atoms. It has recently been shown that in the absence of spin-orbit coupling but 
within an artificial magnetic field a two-component Bose gas forms spin singlets [121 [32] • 
The strongly correlated regime is then described by a series of states with non-Abelian 
anyon excitations (non-Abelian spin singlet or NASS states) [HI H5j. It has also been 
shown that via a Rashba-type spin-orbit coupling a certain spin polarization is favored, 
resulting in spin polarized quantum Hall phases [3H1 H7] which are derived from the 
Read-Rezayi series [35]. Distinct quantum Hall phases, however, have been predicted for 
particular spin-orbit coupling strengths which correspond to a degeneracy between two 
Landau levels [361 [31] • Numerical studies restricted to a single Landau level have been 
able to back the occurrence of the Read-Rezayi-like phases [4*91 150]. For the interesting 
configuration at the degeneracy points, gapped phases at the NASS filling factors have 
been numerically proven recently [42]. 

We undergo a complete numerical study, complementing the results presented 
in [32], thus taking into account two Landau levels. Performing exact diagonalization 
with such an increased basis, our study remains valid also at the degeneracy points. We 
start with a detailed description of the system in Sec. [21 Next, we solve the single-particle 
problem in Sec. El Afterward, we turn to the many-body problem in Sec. HI calculating 
the interaction matrix elements. With this we can perform the exact diagonalization, 
yielding the results presented in Sees. [5] and El Our main findings are summarized in 
Seed 
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2. The system 

We consider a pseudospin-1/2 Bose gas in two dimensions, described by a single-particle 
Hamiltonian 



with I2 being the 2x2 identity matrix. An artificial gauge potential A mimics a 
constant, perpendicular magnetic field and a coupling between the pseudospin and the 
external motion. For convenience, we choose this coupling to be of the Rashba-type. In 
the Landau gauge, the vector potential reads 



Here, B is the magnetic field strength, q controls the strength of the spin-orbit coupling, 
and o~i are Pauli matrices. Both contributions to this gauge potential have already 
been engineered experimentally separately: An artificial magnetic field of this form (i.e. 
q — 0) has been realized through an atom-laser dressing in Bose-Einstein condensates 
without pseudospin degree of freedom [19], and through rotation in two-component 
systems [TTj [T2] . A similar spin-orbit coupling has been achieved experimentally by a 
Raman coupling of three atomic states with two lasers [271 EE]- A theory proposal of 
how to achieve a gauge potential combining both the magnetic and the spin-orbit part 
has been given in Ref. [3B1HI]- Here, the spin-orbit term arises due to a tripod coupling 
of three atomic ground states to one excited state by three laser fields. A review of 
different proposals is given in Ref. |15j . 

All schemes have in common that they effectively consider dressed atoms, that is 
superpositions of different bare atomic states. It is in the basis of dressed states that 
the single-particle Hamiltonian takes the desired form of Eq. ([[]). While equivalent on 
the single-particle level, the different proposals might require different descriptions for 
the particles' interactions. Tripod schemes with different laser configurations, though 
mimicking the same spin-orbit coupling, have been shown to support different mean-field 
phases due to differences in the interaction term [51]. It has also experimentally been 
shown that the interactions between dressed atoms might become long-range, though 
the bare atoms interact via s-wave scattering [28]. Also if the interaction potential of 
the dressed atoms effectively remains a contact potential, the scattering lengths may 
be modified. In the first experimental realization of spin-orbit coupling, the strength of 
interactions between different pseudospin states behaved linearly with the effective spin- 
orbit coupling strength, which caused a quantum phase transition between a pseudospin- 
mixed phase and a pseudospin-separated phase upon tuning the spin-orbit coupling [27] . 

Thus, the choice of interactions is a delicate issue. In this paper, we consider 
two-body contact interactions between the dressed atoms. 



Such a description is appropriate for the pioneering experiment in Ref. [27J. In this 
experiment, 87 Rb has been used with nearly equal scattering lengths in and between 




(1) 




(2) 




(3) 
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different hyperfirie states. Also, for weak or moderate spin-orbit coupling strength 
the scheme of Ref. [27] does not significantly alter the effective scattering lengths 
for the dressed atoms. It is thus convenient to assume SU(2) symmetry, that is, 
fl'tt = 9u = 9n = 9- 

Previous studies have performed numerical diagonalization of pseudospin-1/2 
systems without spin-orbit coupling (q = 0) [121 03]. We note that in this case, also 
the pseudospin polarization S defined as the difference — of up- and down-spin 
particles is a conserved quantity, further reducing the size of the Hilbert space. For 
SU(2)-symmetric interactions, the ground state at almost all filling factors is found to 
be a pseudospin singlet, 5 = 0. For atoms interacting via a two-body contact potential, 
eigenstates with zero interaction energy exist up to a filling factor v = 2/3, where the 
(221) Halperin state becomes the unique zero interaction energy state: Its expression 
on a disk reads [52] 

= n^t - ^-t) 2 nw - zii? n^t - ^ e , (4) 

i<j k<l i,k 

where we have introduced the complex variables, zj = Xj + iyj. 

This Halperin state is part of a series of states occurring at filling factors v = 
2k/ 3, defined as the maximum-filled zero-energy eigenstates of a (k + l)-body contact 
interaction. This series is commonly known as the NASS series, since the states with 
k > 2 support non-Abelian anyonic excitations [H1H5]. As pointed out in Refs. [4*2j |4~3] . 
the non-Abelian phase at filling v = 4/3 is realized even in systems with a two-body 
contact potential. 

Also systems with finite q have been considered numerically [4~9l [50] . The lowest 
Landau level approximation which has been made in these studies is valid within broad 
ranges of q with sufficiently large Landau level gaps. These studies have shown that 
in these regimes the spin-polarized fractional quantum Hall states like the Laughlin 
state [1] or Moore- Read state [3] are relevant, if one maps them into the Landau level 
structure of the spin-orbit coupled system according to Eq. (TI~8l) . We will demonstrate 
this in more detail in Secs.[5]and[6]for the bosonic Laughlin state at filling factor v = 1/2, 
investigating its fate when passing through the degeneracy point q 2 = 3B, and for other 
relevant filling factors which may present incompressible phases, respectively. 

3. Single particle Hamiltonian: Landau level structure 

Prior to solving the many-body problem by numerical diagonalization, it is useful to 
solve the single-particle Hamiltonian. For the spin-orbit coupling of Eq. (j2J), this has 
been done by M. Burrello and A. Trombettoni, yielding a Landau level structure which 
is simply a combination of two Landau levels from the quantum Hall system without 
spin-orbit coupling [4~6l l4"Tj . While they consider a system on a disk, for our numerical 
study it is advantageous to work in an edge-less geometry, as it allows to study bulk 
properties even for small system sizes. We thus choose a rectangle of sizes a and b with 
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periodic boundary conditions, that is a torus, and perform the analog derivation for the 
Landau level structure. 

For convenience, we choose units such that h = 1 and M = 1/2, and define the 
following ladder operators: 

1 



c 



t 



f 2B 
1 



{Px - iPy - iBx), (5) 
■-(Px + ip y + iBx), (6) 



V2B 

with [c, c'] = 1 and [c, c] = 0. The single-particle Hamiltonian can now be expressed in 
terms of these operators 

( 2B{c\ Cl + \) + 2q 2 2V2B gCi \ 

1 \ 2V2Bqc\ 2B{c\ Ci + §) + 2g 2 ) ' Uj 

The diagonal part has the form of a harmonic oscillator, yielding the standard Landau 
levels ^ n)K (x,y) = exp(iny) $„ A (x) and 

$n, K (x) = e- B ^ +K ' B ^' 2 H n [VB(x + k/B)] (8) 

with if n the Hermite polynomials. The energy quantum number n denotes the Landau 
level. Applying the periodic boundary conditions, as done in Ref. [53], k is restricted to 
integer multiples of 2ir/b. This yields a second quantum number j e Z via the definition 
Xj = k/B — 2nj/(bB). The quantity Xj can be interpreted as a displacement of the 
harmonic oscillator in the ^-direction. To ensure the periodicity of the wave function 
in x, with period a, we can sum over all displacements Xj + ka with k an integer, 
V n ,j(x,y) = Y^kL-oc^n,B(Xj+ka)(x,y)- This will not affect the 6-periodicity in y, if 
ah = 2n\ 2 N<$, with A^^ an integer. To prevent double-counting, we must restrict j to 
1 < j < We find that quantizes the magnetic field, and can thus be interpreted 
the number of magnetic fluxes within the unit cell. It is connected to the particle number 
iV through the filling factor 
N p 

where we also have introduced the co-prime integers p and q. Introducing the unitless 
variables x' = x/X and y' = y/X, with A = 1/ y/~B, we finally arrive at the wave functions 
of Ref. [53]. In the lowest Landau level (i.e. for n = 0), they read 



* 0li = (27r 2 iV|-) e-V-^w e 3 



77T / 7T 

J • I ..* Ng,a/b 



(10) 



^ y 2N^a/b 

Here, 6*3 is an elliptic 0-function, and = a; — The excited states can easily be 
obtained by applying the ladder operator c\ \Jn + l\l/ n+ ij = c^\l/ n j. 

Without spin-orbit coupling, i.e. for g = in Eq. (jTJ), the wave function 

*E° = *cu(a|t> + )9|4». (11) 

with |a| 2 + |/3 1 2 = 1 describes the lowest Landau level. Thus, apart from the degeneracy 
in j, we have an SU(2) degeneracy due to the pseudospin-1/2 degree of freedom. 
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Figure 1. Left: First eigenvalues of the single particle Hamiltonian, as a function 
of q 2 /B, see Eq. (|TU|) . Right: Weights of the two components of each eigenvector of 
the single particle Hamiltonian for the "-" states, a~, and (3~ as a function of q 2 /B, 
see Eq. (CG 



For finite g, we have to take into account higher Landau levels, as the off-diagonal 
elements in Eq. ([7]) provide a Jaynes-Cummings-type coupling. Thus, the eigenspinors 
must be of the form 

*^. = ( a ±^- ) /3±^ +1J ) T ) (12) 

where and (3^ are functions of q. The single-particle spectrum is found to be the 
same as in Refs. jHU H7], as expected due to gauge invariance. It reads 

E± = 2B(n + 1) + 2q 2 ± ^ B 2 + 8Bq 2 (n + 1) . (13) 

In Fig. []]we present the lowest eigenvalues E^. 

Since the eigensolutions corresponding to E£ are higher than the lowest E~ 
solutions, we can restrict ourselves to the latter. Their coefficients a~ and (3~ read 

a~ = Af (B + 2q^/2B(n + 1) + \J B 2 + 8Bq 2 (n + 1)) , (14) 
0~ —Af (B - 2 qv /2B(n + l) - y/B 2 + 8Bq 2 (n + l)^ , (15) 
Af = - (b 2 + 85g 2 (n + 1) + 2 gv /25(n + l) 

x ^B 2 + 8Bg 2 (n + 1)) V2 . (16) 

In Fig. [T] we plot the value of the coefficients a~ and f3~ as a function of the non-abelian 
strength, q. The lower spin component quickly acquires weight, measured by |/3| 2 , as 
the spin-orbit coupling is turned on. 

We shall distinguish two cases, as seen in Fig. HJ or Eq. ( fT3l) : 1) For broad ranges 
of q, the ground state is given by a single E~ level. Thus, the pseudospin degeneracy 
of the system with q = is lifted, and the system is effectively spin-polarized. 2) The 
second case is restricted to separate points at q 2 /B = 2n + 3, where two solutions, E~ 
and E~ +1 , become degenerate. This scenario, with a two- fold degenerate single-particle 
ground state, is formally similar to the spin-degenerate case of Eq. ( Illjl . We can define 
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a pseudospin-1/2 spinor as, 

= 7*oj + **rj with | 7 | 2 + |5| 2 = 1. (17) 

The connection between these two Landau levels of the spin-orbit coupled system and 
the q = 0, spin-degenerate lowest Landau levels, is made through the bijective mapping 

|0,j,t) H- * -, = a " \0JA) + /9 " (18) 

|0,jU) ^ *ij = «r |1, J, t> + ft" |2,j,|) • (19) 

For brevity, in the following we denote these two levels by \n = 0, j) and \n = In 
our numerical calculations, we will explicitly take into account these two levels, as we 
focus on configurations where q 2 /B is close to the first critical value of 3. 

4. Many-body Hamiltonian 

From the single-particle states, we easily construct many-body Fock states 

\n n= oj = i, ■ ■ • , n n= o,j=N<s,, n n= ij=i, • ■ ■ , n n =ij=jv$) , (20) 

and we can express the Hamiltonian in terms of annihilation and creation operators, b^- 
and b n j. These obey bosonic commutation relations and act on the occupation number 
n n j in the usual way. The Hamiltonian reads in this notation: 

K = Yl E ™ ^nfinj + Yl V {njAjAijfin 3 j 3 K Aj4 . (21) 
n J {n,j} 

The interaction matrix elements V{ n j}, 

V {n,j} = (ni,h\ (n 2 ,32\V\n 3 ,j 3 ) \n 4 ,j 4 ) , (22) 
with V given in Eq. ([3]), and {n,j} denoting the set of quantum numbers n±, . . . , n 4 and 



ji, . . . , j 4 , are given in Appendix A 



Before numerically diagonalizing the Hamiltonian % of Eq. (T2TT) . we shall exploit 
the translational symmetry of the system on a torus. This symmetry, fully discovered 
by F. D. M. Haldane p3], yields a many-body basis where states are characterized 
by a "Haldane momentum" K = (K x , K y ) which is conserved by the Hamiltonian. 
A comprehensive recipe for constructing this basis can be found in Ref. [55] • The 
main idea is to divide the many-body states into equivalence classes: An A^-particle 
state characterized by the quantum numbers {ji, . . . ,Jn} (and the n quantum numbers 
which here do not matter) can by translation along x be transformed into states with 
{ji — qm, ■ ■ ■ iJn ~ qm}, where q is defined in Eq. (Q, and m is an integer. It runs from 
to some maximum value c v , where . . . , jjv} = {ji — c v m, • • • ,Jn — c v m} due to 
the definition of j modulo N$>. This set of c v different states forms an equivalence class 
which we label by n. We denote each state by \n, m). For each class 77, the translational 
symmetric eigenstates are found as superpositions which can be characterized a quantum 
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numbers J running from to N$ /q — 1 : 

VAIJ)) - ^X>P (—) (23) 
V c ^ fc =o V r i J 

A second quantum number J is given by J = jjmodN^. Both quantum numbers J and 
J can be related to the Haldane pseudomomentum K: 

= -*•*<■'-■«)■ < 24 > 

The quantum numbers of zero pseudomomentum are defined as the point in the Brillouin 
zone with highest symmetry. In most cases, (Jo, Jo) = (0,0), but if N is even and 
p and q are both odd, other values become possible (cf. Ref. [55]). The numerical 
diagonalization is then performed within a Hilbert space with fixed filling factor v and 
fixed pseudomomentum K. 



5. The fate of the Laughlin state at v = 1/2 

The Laughlin state, which on a disk is given by a wave function 

*L=n^-^) 9e_E<i2i|2/4 ' ( 2s ) 

i<j 

is the unique zero-energy eigenstate of a spinless or spin-polarized system of particles 
interacting via a two-body contact potential at v = 1/q. We can thus easily find its 
second-quantized expression by exact diagonalization. Performing the mapping into the 
spin-orbit coupled Landau levels according to Eqs. ( |T8l) and (fT9l) is then a simple matter 
of redefining the basis of the Hilbert space. Since the Laughlin state is a spinless state, 
there is an arbitrariness in applying the mapping of either Eq. ( fl8l) or Eq. ( fT9l . With the 
first choice, we obtain a spinor wavefunction, which in its upper component is exactly 
the Laughlin state, while its lower component is a Laughlin state where all n-quantum 
numbers have been shifted from to 1. With the second choice, this shifted Laughlin 
state describes the upper component, and the lower one is a Laughlin which has been 
shifted into the second Landau level. 

Let us investigate in which regimes of the spin-orbit coupling strength such 
projected Laughlin states describe well the ground state at filling factor v = 1/2. 

5.1. Below the first degeneracy point: < q 2 /B < 3 

For q 2 < 3B, the levels \n = 0,j) are the lowest single-particle states, see Fig. [TJ 
thus we expect the ground state to be given by the projection into this manifold, 
following Eq. (|T8|) . As shown in Fig. [2] (c), the overlap of the true ground state with 
this projected Laughlin state is approximately 1, if we are sufficiently far below the 
degeneracy point. However, it quickly decreases to when approaching the degeneracy 
point. Simultaneously, the n = Landau level is depleted, see Fig. [2] (b). 
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Figure 2. Properties of the system as a function of the strength of the spin-orbit 
interaction, q 2 /B. (a) Low-lying excitation energies, (b) (left scale) Interaction energy 
of the ground state, (right scale) Occupation of n = level in the ground state, (c) 
Overlap of ground state with different test states as further explained in the text: 
Laughlin projections, and ground states of two-body contact interaction in the first 
and second excited LL. The system size is N — 6 at half filling, with spin-independent 
interactions g SlS2 = 0.2. Energies are in units of gB. 



In Fig. [2] (a), we plot the low-lying excitation energies as a function of q 2 j B. This 
shows that the projected Laughlin phase is gapped. The gap, however, vanishes when 
approaching the degeneracy point. In Fig. [2] (b), the interaction energy of the ground 
state is plotted. The projected Laughlin state has a finite interaction energy, due 
to interactions within the lower component and between the components. A sudden 
decrease in the interaction energy is observed when the degeneracy point is approached. 

The numerical results shown in Fig. |2] reflect the situation of SU(2) symmetric 
interactions. The value for the interaction strengths defines the interval in q where the 
system is occupying the n = 1 level. For the relatively small value g = 0.2 in Fig. [2] (c), 
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the drop in the overlap curve appears to be quite sharp. We can extend it over a larger 
region by choosing a stronger interaction. We have also investigated an interaction 
term, where g^ = g^ and g^ = 0. Qualitatively, this yields the same results as in the 
SU(2)-symmetric case for < q 2 < 3B. 

5.2. At the degeneracy point and slightly above: 3 < q 2 / B < 3.2 

Slightly above the degeneracy point (see Fig. CD), when most particles populate the n = 1 
level, the projection of the Laughlin state according to Eq. ffl9|) represents reasonably 
well the true ground state. As seen in Fig. [2] (a), also this projection yields a gapped 
state. As in the standard Laughlin state, the ground state is unique (apart from the 
center-of-mass degeneracy), and characterized by a Haldane momentum K = 0. This 
gives another hint for the Laughlin-like behavior of the system. 

The projected Laughlin states, however, fail to describe the system precisely at the 
degeneracy point. Here, the gap vanishes and a quasi-degenerate manifold of several 
states at different K describes the ground state. Apart from the fact that K = does 
not belong to this manifold, the lowest state at this pseudomomentum has no overlap 
with the projected Laughlin states. The failure of the Laughlin projections at the 
degeneracy point can be traced back to the fact that the projected Laughlin states are 
restricted to a single Landau level. We expect that at the degeneracy point, the particles 
have to be mapped partly according to Eq. ({TBI , and partly according to Eq. ( fl9|) . but 
we have no unique prescription for this. At the degeneracy point, our numerical results 
are very sensitive to particle number with respect to the number and Haldane momenta 
of the states within the ground state manifold. 

5.3. Clearly above the degeneracy point: 3.2 < q 2 /B < 4 

The regime where the n = 1 Laughlin state projection describes well the system turns 
out to be restricted to a small parameter range above the degeneracy point. When 
further increasing the spin-orbit coupling strength, the overlap between true ground 
state and projected Laughlin state rapidly decreases and finally vanishes. For our 
choice of interaction of Fig. [2] (c), SU(2) symmetric with g = 0.2, it is zero already 
for q 2 /B = 3.2. Here, the system is still far away from the next degeneracy point at 
q 2 /B = 5, thus the two lowest Landau levels are still given by n = 1 and n = 0. 
Taking into account these two levels, we have performed the exact diagonalization for 
up to iV = 6 particles. Since in that region, roughly located around 3.2 < q 2 /B < 4, 
a significant gap separates the n = 1 Landau level from the n = level, we also have 
performed exact diagonalization studies neglecting the n = level, which indeed has 
yielded the same results. With this reduction of the Hilbert space, we have reached 
system sizes of up to iV = 8 particles. 

Again we have considered both an SU(2)-symmetric interaction and the 
configuration with g^ = g^ and = 0. Qualitative differences in the spectral structure 



FQH states of a Bose gas with spin-orbit coupling 



12 



N=6 



— 

c 

o 

■ i— ( 



0.03 
0.02 
0.01 


0.03 
0.02 
0.01 


0.03 
0.02 
0.01 








4+ ,, #4 



+ 44- 4+ +- 



i # — i r 

44- 



"t^ Jit X^-S 



# 



j i i i_ 



3 
1X1 



Figure 3. Excitation energies of low-lying states at q 2 / B = 3.6 for SU(2)-symmetric 
interactions (g = 0.2). Energies are in units of gB, and \K\ is in units of 2irh/a. 



are found for N < 5, but they are absent in larger systems. We thus restrict ourselves 
to presenting the data obtained for the SU(2)-symmetric case. 

For all particle numbers, a state at K = is the ground state of the system, but 
this state is not well described by the projected Laughlin state. Neither is it protected 
by a significant gap: Depending on the particle number, states at different K are 
energetically very close to the ground state, see Fig. |3j The upper spin component of 
this state is found to be well described by the ground state of a one-component system 
which is assumed to live in the first excited Landau level, see Fig. [2] (c). This overlap 
further increases: For N = 8 it reaches the value 0.96, if we renormalize the upper 
spin component to one. Another state which has a sizable overlap with the true ground 
state is the ground state of a one-component system in the second excited Landau level. 
Comparing this state with the re-normalized lower spin component of the true ground 
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state, the overlap is 0.38 for N = 6 and decreases to 0.30 for N = 8. 

The relevance of these states is not very surprising, given the fact that the spin-orbit 
coupled level \n = l,j) is constructed from the first and second excited Landau level 
of the spinless problem. But it is noteworthy that the projection of the Laughlin state 
into these levels has basically no overlap with the respective ground states. Studying 
the behavior of particles in higher Landau levels is particularly relevant with respect to 
fermionic systems. In the field of strongly correlated electrons, very much attention has 
been paid to the Hall plateau at v — 5/2. For a spin-polarized system, this represents 
the situation of a half-filled first excited Landau level, which might be described by 
a Moore-Read-like wave function [3] and thus support non-Abelian excitations. With 
bosonic systems, it is not possible to study higher Landau levels by filling the lower 
levels. However, as our results clearly prove, the spin-orbit coupled system above the 
degeneracy point behaves with very high fidelity like a one-component system in the 
first excited Landau level. It thus becomes possible to explore this regime also within 
bosonic systems. 

From the results shown in Fig. [31 it is unclear whether the system develops a gapped 
phase in the thermodynamic limit, as the tendency observed for sizes up to N = 7 may 
suggest. At N = 8, however, a quasi-degeneracy of the K = state with several states at 
distinct Haldane momenta, indicates that a phase with a broken translational symmetry 
is established, like crystalline phases or bubble phases (56J [571 [58] - Such phases have 
been discussed in the context of fermions as possible candidates for substituting the 
Laughlin state in higher Landau levels. 

6. Incompressible phases at the first degeneracy point 

We have seen in the previous section that the system behaves as an incompressible 
liquid at v — 1/2 for a spin-orbit coupling strength extending from < Q 2 /B < 3, 
and 3 < q 2 / B < 3.1. These phases could be understood as Laughlin phases projected 
into a Landau level structure which takes into account the spin-orbit coupling. Such a 
mapping fails at the degeneracy point, q 2 /B = 3, as the number of states within the 
lowest spin-orbit coupled Landau level is doubled. Accordingly, the incompressibility 
for v = 1/2 is lost. In this section we investigate the system at this degeneracy point, 
and look for filling factors which might support gapped phases. 

6.1. u = 2 

A clearly gapped phase shows up at filling v — 2, where a unique ground state at K = 
is found for N = 8 and N = 10. This fact makes this phase distinct from the Read- 
Rezayi states and the NASS states which are both possible candidates due to the filling 
factor. Apart from the large gap, the robustness of the spectrum against deformations 
of the torus ratio, shown in Fig. HJ signals the incompressibility of this phase. 
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Figure 4. Spectrum at v = 2 for SU(2) symmetric interactions. Axis ratios 1 and 
0.75 are plotted. Note that the degeneracy between 90° rotations in it-space is lifted 
for ratios different from 1. and therefore, more states seem to be present at ratio 0.75. 
Energies in units of gB, and \K\ in units of 2nh/a. 



6.2. v = 3/2 

As has been pointed out in several works [101 091 [50] , off the degeneracy points 
incompressible phases are supported not only at v = 1/2, but also at other fillings 
corresponding to the Read-Rezayi series |48j. At the Read-Rezayi filling, v = 3/2 for 
N = 6, a ground-state at K — is separated by a sizable gap of about twice the 
typical energy difference from a second state at K = 0. Since a two-fold degenerate 
ground state at K = characterizes the v — 3/2 Read-Rezayi state on the torus, one 
might expect that at this filling a Read-Rezayi phase exists on the degeneracy point. 
We have thus calculated the ground-states of a spin-polarized system with four-body 
contact interaction. Projecting the two zero-energy states into the LL structure of the 
SO coupled system, we obtain an overlap of 0.70 with the real ground state. This 
number is close to the total weight of the state within the polarized subspace P = 6, 
being 0.72. For the second state, an overlap of 0.32 is found. Increasing the particle 
number to N = 9, we still find overlaps of 0.39 and 0.42, but the spectral structure is 
not robust. While one state at K = remains the ground state, the second state at 
this pseudomomentum lies above other states in the spectrum. 

6.3. v = 1 

In spin-polarized Bose systems, the famous Moore-Read state is supported at filling 
factor v = 1. Also for spin-orbit coupled systems, this state plays a role: Below the 
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Figure 5. Spectra of a two-component Bose gas without spin-orbit coupling (q = 0) 
at filling factor v = 2/3 for different spin polarizations S. N = 6. 



degeneracy point (q 2 < 3B), the three lowest eigenstates are found at K = (iV$/2,0), 
(0,iV$/2), and (iV$/2, iV$/2). The phase can thus be identified with the Moore- Read 
phase jl9]. However, this spectral structure disappears at the degeneracy point. Here 
we find ground states at different pseudomomenta, which even depend on the size of 
system: (3,3), (3,5), (5,3), (5,5), (4,0), and (0,4) for iV = 5, and (3,0), (0,3), (2,2), 
(2,4), (4,2), and (4,4) for N = 6. For N = 7, states at (2,0), (5,0), (0,2), and (0,5) 
provide the ground state manifold. For N = 8, we have (3, 3), (3, 5), (5, 3), (5, 5), (4, 0), 
and (0,4). Since the gap between these states, compared to other energy differences in 
the spectrum, is not extraordinarily big, the system is expected to be compressible. 

6.4- NASS series 

The NASS filling factors, v = 2/3 and v = 4/3, have been recently shown to support 
incompressible phases in the system without spin-orbit coupling [42} 03] . It is instructive 
to recall the results in that limit, q = 0, for our subsequent discussion at the first 
degenerate point. In Fig. [5] we depict the spectrum of a system of 6 atoms interacting 
through a two-body contact interaction at q = at filling v = 2/3. While for |5*| = 6 and 
\S\ = a state with K = is the unique ground state, a (quasi) degenerate manifold of 
states with K ^ provides the ground state for spin polarizations \S\ — 2 and \S\ =4. 
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Figure 6. Spectra of a two-component Bose gas without spin-orbit coupling (q = 0) 
at filling factor v = 4/3 for different spin polarizations S. The results for a two-body 
contact interaction are given by the black crosses which correspond to the energy scale 
on the left. The red triangles with the energy scale on the right show the results for a 
system interacting via a three-body contact potential. 



Interestingly, the energy of low-lying excitations at \S\ precisely agree with the ground 
state energy at |S"| = |>S| + 2. 

Similarly, for the is = 4/3 filling, we have compared the spectrum obtained for 
a two-body contact interaction with the one obtained with only three-body contact 
interactions. Both spectra are depicted in Fig. |6j The experimentally feasible case of 
two-body contact interactions shares most of the spectral properties found in the three- 
body contact case in all polarizations. In particular, it is especially telling to note that 
the degeneracies of the lowest states for each polarization are shared in both calculations, 
a signal of the topological equivalence between both systems [I]. For S = 0, the gapped 
ground state of the three-body contact interaction defines a NASS phase. 

Now, we turn back to the spin-orbit coupled system. Here we do find a clear gap 
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Figure 7. Spectrum at the degeneracy point q 2 /B = 3 for filling factor ^ = 2/3 with 
iV = 6 (left) and for v = 4/3 with N = 8 (right). We have chosen SU(2) symmetric 
interactions. Energies in units of gB, and \K\ in units of 2nh/a. 



in the spectrum at v — 2/3 with AT = 6, and a small gap at ^ = 4/3 with N = 8, see 
Fig. [7J 

For the u — 2/3 filling, the nature of the ground state of the system can be related 
to the 221-Halperin state, Eq. fl3j). We have projected the 221-Halperin state, onto 
the Landau level structure of the spin-orbit coupled system according to Eq. (TIB"]) , and 
computed the overlap of this state with the v = 2/3 ground state. We find the value of 
0.18. It should be noted, however, that interactions now can induce "spin" flips within 
the two-fold degenerate manifold. Therefore the analog of spin polarization, namely 
population imbalance P = No — N\, is not conserved. Here, N\ (No) define the number 
of particles in the n — 1 (n — 0) Landau level. We therefore write the true ground state 
|^) as 

JV 

|W)= u>i\P = i) 9 , (26) 

i=—N 

with \P = i)^ the contributions with well-defined population imbalance. Mapping the 
Halperin state with spin polarization S — to the spin-orbit coupled Landau levels 
according to Eq. (TlBl yields a state with P = 0. The overlap of this state with 
|\&) thus is bounded by the corresponding weight, wq = 0.26. Contrarily, the overlap 
of the projected Halperin state with \P = 0)^ takes a significantly increased value, 
0.18/0.26 w 0.7. 

Interestingly, the most important contributions to \^) are \P = 6)^ with a weight 
w 6 = 0.64, \P = 4)^ with a weight = 0.57, and \P = 2)^ with a weight w 2 = 0.42. 
The contributions with P < are rarely populated with weights W-2 = 0.1, W-4 = 0.07, 
and u>_ 6 = 0.01. This shows that population of n = is still favored, despite the 
degeneracy with n = 1 on the single-particle level. It reflects a fact which we have 
already encountered before when analyzing the system at v — 1/2: Contact interactions 
favor the Landau level with lower n. 

As mentioned above, with good overlap we can interpret the P = contribution as a 
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projected Halperin state, that is, as the ground state of a two-component system without 
spin-orbit coupling (q = 0) in a two-body contact potential with spin polarization S = 
and filling factor v = 2/3. Diagonalizing such a system for different S, and projecting 
the lowest state with K = according to Eq. (fl8|) . we are able to relate all the states 
\P)q in the decomposition of Eq. (1261) to ground states of the two-body contact potential 
for a fixed spin polarization S and Haldane momentum K . Explicitly, we find an overlap 
of 0.94 between the projected ground state at S — 6 and \P = 6)^, an overlap of 0.84 
for S = 4, and an overlap of 0.83 for S = 2. Modeling a state with the projected ground 
states for all spin polarizations — 6 < S < 6, we are able to reproduce the true ground 
state with a fidelity of 0.82. 

For v = 4/3 with N = 8, the overlap with the projected NASS state is only 0.03. 
Nevertheless, the true ground state can be reproduced from the lowest K = states of 
a three-body contact interaction, if we take into account different spin polarizations S 
and project them according to Eq. ( ITS"]) . With this we are able to model the true ground 
state of the spin-orbit coupled system with a total fidelity of 0.52. 

6.5. v = 2/9 

Finally we note that an incompressible phase can be expected at critical uq, defined as 
the largest filling supporting states with zero interaction energy. For any v < uq, the 
interaction energy should remain zero, while for v > Vq it is by definition larger than 
zero. Thus, at z/o, a kink in the energy as a function of v can be expected, resulting in 
an incompressibility of the system. As Fig. [2] (b) shows, the spin-orbit coupled system 
with two-body contact interaction has a finite interaction energy at filling v = 1/2, in 
contrast to the system with only an Abelian gauge field. This is a consequence of the 
fact that, at this filling, there are no zero-energy states of two-body contact interaction 
within the first excited or higher Landau levels. We thus should focus on smaller filling 
factors. 

By analytical arguments, a zero-energy state has been predicted at v = 2/9 in 
Ref. [J7]. Such small filling factors are difficult to reach numerically by means of exact 
diagonalization, as even a small particle number requires a relatively large number of 
fluxes. We have performed the numerics for N = 4 and iV$ = 18, and found, for each 
of the allowed pseudomomenta (0,0), (1,0), (0,1), and (1,1), exactly one state with 
negligibly small interaction energy. In contrast to this, the exact diagonalization at 
v = 1/4 (with N = 4) yields several zero-energy states, whereas for v = 1/3 (with 
N = 4 and N = 5) no zero-energy states can be found anymore. 

However, we should also note that the energy as a function of v is concave. In the 
thermodynamic limit, it might happen that ^ — > 0, when v \ |. This would render a 
compressible phase even at this filling. We therefore conclude that, at the degeneracy 
point, robust incompressible phases are rather located in the denser regime. This should 
facilitate the experimental observation. 
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7. Summary and conclusions 

We have performed a numerical study of the ground and excited states of an 
ultracold pseudospin 1/2 bosonic cloud subjected to an artificial magnetic field. The 
setup considered, first proposed in Ref. (46], includes a spin-orbit coupling, which is 
mathematically equivalent to a non-abelian SU(2) gauge potential. The single particle 
spectrum has a number of degenerate points which have received a special attention, 
as they mix different Landau levels of the SU(2) symmetric Hamiltonian, a fact which 
makes them particularly appealing for the existence of unknown interesting new phases. 

We have concentrated on the parameter region around the first degenerate point 
q 2 /B = 3. First, we have analyzed v = 1/2 and discussed the different Laughlin-like 
states which have a large overlap with the ground state of the system. The exact ground 
state of the system at this filling is described mostly either by the usual Laughlin wave 
function (mapped into our Landau level structure), if we are below the degenerate point, 
or by the corresponding ground state of contact interactions in the first excited Landau 
level or second excited Landau level, for values of q above the degenerate point. This 
allows to explore higher Landau levels even with bosons. Slightly above the degenerate 
point there is a sizable overlap between the ground state and a Laughlin state projected 
into the first excited Landau level. 

Secondly, we have studied in detail the phases appearing at the first degenerate 
point for the most interesting filling factors. Notably, we find an unexpected gapped 
phase at K = at a fairly dense filling, v — 2. This phase is gapped for N = 8 and 
N = 10 and is robust for changes in the aspect ratio of the torus, as expected for an 
incompressible phase. Motivated by the v — 1/2 results, we looked for incompressible 
phases at other Read-Rezayi fillings, v = 1 and v — 3/2. For v = 1 we find no indication 
of a gap in the spectrum of the system. For the v — 3/2 we find a sizable overlap between 
the exact ground state of the system and the ground state obtained with four-body 
contact interactions mapped into our Landau level structure. The spectral structure 
however is not robust as we increase the particle number, finding that the two-fold 
degeneracy of the K = ground state is lifted when we go from 6 to 9 particles. 

Thirdly, we have analyzed the gapped phases occuring at the NASS filling factors 
v = 2/3 and v = 4/3. Decomposing the ground states into the contributions with a fixed 
population imbalance between the two Landau levels, we are able to relate these states 
to the ground states of a two-body contact potential for v = 2/3, or of a three-body 
contact potential for v = 4/3, obtained for fixed spin polarizations in a system without 
spin-orbit coupling. 

As a conlusion, we note that the properties of the spin-orbit coupled system in 
the vicinity of the degeneracy point make the considered setup very appealing from the 
experimental point of view: Upon tuning the spin-orbit coupling strength, transitions 
between different topological phases can be realized. It also becomes possible to explore 
the behavior of bosons in higher Landau levels. For the relatively dense system at filling 
v = 2, a robust gapped phase occurs which is neither described by the NASS nor the 
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Appendix A. Interaction matrix elements 



The interaction matrix elements Vf n> j\, with {n,j} denoting the set of quantum numbers 
rii, . . . , n 4 and ji, . . . , j 4 , are given by 

V {n,j} = ( n uh\ {™2, 3z\V \n s , j 3 ) \n 4 ,j 3 ) , (A.l) 
with V given in Eq. ([3]). This yields: 

Vs \ — -( a-M- a~ a~ a~ a~ X r - J • 



r»4 



+ 9U fin$ni$nzfin A ^nf* 



+I,n 2 +l,n 3 +l,n4+l 



rji J2 ,33 J4 



"I" 9ti a n 1 Pn2 a n 3 (^n4 -^ni,n 2 +l,n3,n4+l 

I „ O— — O- — T-jl ,J2j'3J4 

y^t l J n 1 Lx n 2 l J n 3 Lx n4 - t -ni+l,n2,n 3 +l,n4 /' 



(A.2) 



with 



X, 



{/} 



{»} 



dVidV 2 5{ri-r 2 )y* nijl (ri)y* n2j2 {r 2 ) 

Jo 

X ^3.73(^)^4.74(^2). 



(A.3) 

To evaluate this integral, we write the elliptic 9 functions in ty n j as infinite sums, see 
also [551 [59]. We get 



-{n} 



2ttN, 



it Yl Yl 5 'h 



+32 ,33+34 jl-ji,t 



fi=—oo u=—oo 



X C, 



x exp 



-r 



) 2 + M 2 ) +irth-h) 



where 5' is a Kronecker delta modulo r = \/lix /N$, and $ = -^/a/ 
Cn t ,nj read (with «,nGK: 

Ci |0 (m,i;) = - (u + iv)/V2= -C% tl (u,v), 
C 1A (u,v) = 1 + C ls0 (u,v)C 0tl (u,v), 
C 2fl {u,v) = V2[C 0A (u,v)} 2 = -Cl >2 {u,v), 



(A.4) 

The coefficients 
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Ci, 2 {u,v) = (u 2 + v 2 -±)(u + iv)/± = -C£ tl (u,v), 
C 2 ,2(u, v) = (8 - 8(u 2 + v 2 ) + 2u 2 v 2 + u 4 + v 4 )/8. 

The sums in Eq. (1A.4|) converge quickly, especially within the lowest Landau level. For 
levels up to n = 2, taking into account only values with < 16 and \u\ < 8 provides 
very precise values. 
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